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Abstract: Recently, Giraitis et al. (2003, [10]) proposed the V/S statistic 
for testing long memory in random sequences. Wc generalize this statistic 
to the setting of random fields. The null hypothesis is concerned with short 
memory random fields while the alternative contains a very large family of 
long memory random fields. Contrary to most of the previous works dealing 
with long-range dependence, no assumption is made about the isotropy of 
the strong dependence. Some simulations are presented in order to assess 
the power of the test according to the kind of long memory in presence. 
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1. Introduction 

A stationary random field X = (X„)„g^£i is usually said to exhibit long mem- 
ory when its covariance function n S Z'', is not absolutely summablc: 
SnGZ"* kl"-)! ~ alternative definition relates on spectral properties: A 
random field is said to exhibit long memory if its spectral density is unbounded. 
These two points of view are closely related but not equivalent. In this paper, 
the concepts of "strong dependence" and "long-range dependence" are the same 
as "long memory" . 

Most of the previous studies on long memory random fields (see [G, 7, 18]) 
assume that the strong dependence occurs with the same intensity in all direc- 
tions. Indeed, these works are concerned with isotropic long memory according 
to the following definition. 

Definition 1. A stationary random field exhibits isotropic long memory if it 
admits a spectral density which is continuous everywhere except at where 

where ||.|| denotes the Euclidean norm, where L is slowly varying at infinity and 
6 is a continuous function on the unit sphere in . 
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However, it is easy to construct non-isotropic long memory random fields. 
In [14], such fields arise from particular filterings of a white noise, from the 
aggregation of weakly dependent random fields, or from systems of statistical 
mechanics in phase transition. 

Our aim is to construct a procedure for discriminating between weak depen- 
dent random fields and strong dependent ones, regardless of the isotropy. 

In dimension d ~ 1, several tests for long memory are available. They are 
mostly based on an estimation of the variations of the partial sums process of 
X. For all these tests, the alternative hypothesis consists in parametric families 
of long memory processes, typically FARIMA time series. Lo (1991, [17]) first 
developed a test based on the R/S statistic, which estimates the range of the 
partial sums process of X. The KPSS test was initially developed by Kwiatkovski 
et al. (1992, [12]) for testing stationarity (under weak dependence assumptions) 
against the presence of a trend or a unit root. A variant of the KPSS test, based 
on an estimation of the second order moments of the partial sums of X, was 
proposed by Lee and Schmidt (1996, [IG]) in order to test long memory. From 
the same idea, Giraitis et al. (2003, [10]) introduced the V/S statistic, based 
on an estimation of the variance of the partial sums of X. It appears that this 
test is more powerful than the R/S test and than the KPSS test for detecting 
strong dependence. 

Note that the behavior of the partial sums of X is not the unique way to test 
long memory in dimension d = 1. Goodness-of-fit tests for long range dependent 
time series have been developed (cf. [1, 9] and [4]). To the best of our knowledge, 
in dimension d > 1, there exists no generic model able to exhibit different situa- 
tions such as isotropic or non-isotropic long memory. A goodness-of-fit approach 
seems therefore too restrictive in our framework. 

In this paper, we generalize the V/S test to the setting of random fields. 
Since this test is based on the partial sums oi X, we summarize in Section 2 
what is known about their limiting behavior under short memory and long 
memory. This will allow us to specify our testing hypothesis. Section 3 proves 
the consistency of the test through asymptotic results. Some simulations are 
presented in Section 4 when d = 2. They reveal that the power of the test is 
strongly related to the anisotropy of the long memory. 

2. Test hypotheses 

Let X be a second order stationary real random field. We want to test the 
null hypothesis: X is weakly dependent against the alternative assumption: X 
exhibits long memory. The V/S test that we extend in this paper does not rely 
exactly on these testing hypotheses. It rather focuses on the behavior of the 
partial sums process, defined for all i G [0, 1]'' by 

E 
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with An{t) = Z'^n riiLi [Ij [('^ ^ l)*iJ + l]j where [ccj denotes the integral value 
of X and n is a positive integer. Indeed, there is a close relation between the 
dependency of X and the asymptotic behavior of its partial sums. 

When X is a centered and weakly dependent random field, i.e. when its 
covariance function is absolutely summable, it is well known that if tr^ := 
X]/iGZ<* '''W 7^ ^ functional central limit theorem generally holds: 

' ^ X^^^^'^Bit). (2.1) 



Here, B denotes the Brownian Sheet, i.e. the centered Gaussian process such 
that E{B{t)B{s)) = HiLi ^ ■s*- The convergence takes place in V{[0,lY), 
the Skorokhod space of cadlag functions defined on [0, 1]'' (see [2] for instance). 
This result has been proved by different authors according to the kind of weak 
dependence of X, among others: Wichura (1969, [-'■ ]) when X is an i.i.d process, 
Dedecker (2001, [5]) under a weak projective assumption. 

On the other hand, when X is long-range dependent, (2.1) is generally false. 
More precisely, when the spectral density of X exhibits at least one singularity 
located at zero, then, under some structural hypotheses on X, it has been proved 
that 

J, E X.-<i^HlV(t), (2.2) 

where 7 > d/2 and y is a random field different from the Brownian Sheet 
(not even necessarily Gaussian). This result is shown in [d] for functionals of 
Gaussian fields, extended in [IS] to functionals of linear fields, when the long 
memory is isotropic according to Definition 1. In the case of linear fields which 
exhibit anisotropic long memory, (2.2) is proved in [15]. 
Let us now precise the testing hypotheses. 

HO: Short memory hypothesis. The second order random field X is sta- 
tionary, with a covariance function r, such that SMI, SM2 and SMS below are 
satisfied. 



SMI 



SM2 



^ |r(j)| < 00 and := ^ r(j) > 0. (2.3) 



fceA„(t) 

where B is the Brownian Sheet. 



5: iX,-EiXo))'''^>Bit), 



SM3 



sup ^ \c4(i,j,k)\ < 00, 



where C4 represents the fourth order cumulants of X: Denoting Xi = Xi ^ 
E{X,), Ci{i,j,k) = E[XoX,XjXk]-r{i)r{k^j)-r{j)r{k-i)-r{k)r{j-i). 
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HI: Long memory hypothesis. The second order random field X is sta- 
tionary and satisfies LMl and LM2 below. 



LMl 



1 ^ T^/V \\ ■^(["'1]'') 



E iX,-E{Xo))'''^>Yit), (2.4) 



with 7 > (i/2, where i is a slowly varying function at infinity and Y is 
some measurable and non-degenerated random field. 

LM2 

Varf J2 X^ =Oin^''L^in)). (2.5) 

Remark 1. The null hypothesis HO deals with short memory as suggested by 
assumption SMI. Assumption SM2 claims that a functional central limit the- 
orem holds, as expected in the short memory setting. However, SMI does not 
imply SM2. For instance, some counter-examples are available in Theorems 7 
and 8 in [11]. Finally, SMS is needed for technical reasons. 

As explained above, (2.4) holds for a large family of long memory fields 
and it can be reasonably chosen as the alternative hypothesis. However, when 
the strong dependence involves non zero spectral singularities, then (2.1) may 
remain true (cf. [ I 'l], Theorem 2 and Theorem 4). Then, there is no chance 
for this particular situation of long memory to be detected by the test. Notice 
that the same restriction exists in dimension 1 in all the long memory testing 
procedures based on the behavior of the partial sums. Assumption LM2 is 
convenient for technical reasons and appears to be a weak restriction to LMl. 

3. The testing procedure 
3.1. Test statistic 

We generalize the V/S statistic to d > 1. Let us first introduce some notations. 
For all positive integer n, let An = ^ti(I)- Wc denote, for all positive integer j, 

S*n,,= J2i^^-^r^), (3.1) 

where Xn = n^'^'^j^j^ Xj. Let q be an integer in [l,n]. An estimator of cr^, 
defined by (2.3), is 

sl= E (3-2) 



with Bq = {—q, . . . , q}'^. Here Wqj = ni=i(l ~ q) some weights leading to 
the positivity of (see for instance [3] p360) and f is the empirical covariance 
function: 



^ n-\ji\ n-\jd\ 

^^■^^ = ^7 ■ ■ ■ i^ku-,k„ - Xn) {Xk, + \j,\....,k^ + \j^\ - Xn) 
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The statistic V/S is defined by 



M„ ^ n 



where Var {S^ j , j G A„ 
One can rewrite M„ as 

7 

Mn = - 



-2d 



(3.3) 



3.2. Consistency of the test 

In what foUows, the integer q involved in Definition 3.2 is actually a function 
Qn depending on n. The following proposition establishes the consistency of the 
test when the sequence Qn is properly chosen. 



Proposition 1. // lim„ 
(i) Under no, 



,qn^oo and lim„_,oo ^ti/"- ~ 0; then, 



(3.4) 

denotes the convergence in 
(3.5) 



where B is the Brownian Sheet on [0, 1]'' and — 
law. 

(ii) Under HI, 

p 

Mn * OO, 

p 

where > denotes the convergence in probability. 

The testing procedure is the following: Given a significance level a G [0,1], 
one rejects the null hypothesis HO if M„, given by (3.3), is greater than c(a), 
where c(a) is such that 

P{Ud > c(a)) = a, 

where Ud is distributed according to the asymptotic law involved in (3.4). Propo- 
sition 1 guarantees that the significance level of the test is asymptotically correct 
and that the power of the test goes to 1. 

Remark 2. To our knowledge, the theoretical form of the asymptotic law Ud 
in (3.4) is unknown for d > 1. When d = 1, the distribution function of Ui is 
Fk{t^\/x), where Fk denotes the Kolmogorov distribution function (see [10]). 
This identification comes from Watson (1961, [19]) and is not easy to extend 
to d > 1. It is however straightforward to obtain E{Ud) = (1/2)'' + (1/4)*^ - 
2(1/3)'*. For d = 2, a simulation of the density distribution of U2 is presented 
in Section 4.1 (cf. Figure 1). 
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Proof. Notice that if ^ < ii < then 

^^(LjiiiJ + 1, . . . , lntd\ + 1) = 5,:(fc + 1, [niaj + 1, . . . , lntd\ + 1), 
where S'*(j) = S** ^ is defined for all j £ An in (3.1). Therefore 

Yl J = [ S^ilnhl + 1, . . . , lntd\ + l)dt. 



The same equality holds with respect to S** and expression (3.3) of A'/„ becomes 



M„ 



[OA]'' 



[0.1]-' 



S'*(LntiJ +l,...,lntd\ +l)dt 



Let Snit) = n '^/^ J2k<£A„{t) ^fe' then Af„ can be expressed as 



/ (sn+iit) ~ n 
'if ^n+i(i)-n 



lnti\ + 1 



dt 



[nUl + 1 



2 -I 



Sn{l)dt 



Hence s^M„ is of the form (f>(S'„(.)), where $ is a continuous map. 

As a consequence, under SM2 and from the continuous mapping theorem, 



Besides, under LMl, 



'n27L(n)- 



The proof is concluded thanks to the following lemma. 
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Lemma 1. // lim„^oo Qn 
(i) Under no, 

(a) Under HI, 



oo and lim„^oo Qn/^^ = 0, then, 



-2 P 2 
s„ — > a . 



rsl ^ 0. 



n27L(n)2"" 
Proof of Lemma 1. 

The demonstration is an adaptation from Giraitis et al. (2003, [10]) to the 
random field framework. 

Denote k = {ki, . . . , kd) and \j\ = (|ji|, . . . , \jd\), 



^ n-\ji\ n-\jd\ 

fci=i fed=i 

where fi is the expectation of X. We spHt 

4 = XI '^9J^(-?')+ 5Z ^q,3{^{j)-r{j)):=Un + Vn. (3.6) 

We first show that, under HO, i?(|u„|) — > when n oo, while under HI, 
i<;(|t!„|) = o(iT?"i~'^L{nY) when n — > cx). 
Some computations lead to 



r{j)-r{j)=mi-^-f]{X 



\ II I 

1=1 

X • • • X ((Xfe - m) + [Xk+\,\ - /.)) 
From the Cauchy-Schwartz inequality 

i?(KI)< E i?|r(jV^~0-)l 

ie-B,-i 



fci = l kj = l 



\ 



fel = b'i| kd = \jd\ 



Since r is bounded, if Si < i^, for alH = 1 . . . d, 

^ E---E(^^-/^) = E ••• E K*-o<cn(^.-^o, 
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where c is a positive constant. Therefore 

E{K\)<c ^ (e (X„ - m)' + 2n-'^y's (X„ - m)' n - ) 
ies,_i V i=i / 

<c J2 + 2n-'/^^E{X^ -l^f)- (3.7) 

Under SMI, E (Xn - l^t)'^ < cn-"^, while under LM2, £'(X„-^)^ < 
cn^^~^'^L^{n). Since q/n goes to 0, it is then straightforward to conchide that, 
under HO, vanishes and, under HI, iJdwnl) = o{n'^''^^L'^{n)). 

Now, we turn to the asymptotic behavior of u„. 

Let us first show that under HO, it„, defined by (3.6), converges in probabihty 
to when n goes to infinity. Notice that 

According to (3.6), converges in probabihty to cr^ if E (ii„ — E{un))^ con- 
verges to 0. 

E {un - E{Un)f 

= '^9:i'^9J'C0v(f(j),f(j')) 

< E |cov(r(j),f(j'))| 

We split the sum above in two terms involving on one side the cumulants and 
on the other side the covariance function. Indeed 

cov {{Xk - Ai)(Xfe+|j| - /i), {Xk' - fi){Xk'+\j'\ - Ai)) 
= cMm{Xk, Xk+\j\, X'^, Xk'+\j'\) 

+ r{k - k')r{k' -k+ l/l - IjI) + r{k - k' - \f\)r{k - k' + 
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First, thanks to SMS, 

H H \cum{Xk,Xk+iji,X'^,Xk>+ifi)\ 

- H H H \cum{Xa,Xyi,X„X,,)\ 

H E MIj|,m')I 

where c is a positive constant. 
Then, 

^ E E \r{k-ky{k'-k+\f\-\j\) + r{k-k'-\f\)r{k-k' + \j\)\ 
E E E 2|r«r(z')| 




where c is a positive constant. 

Finally E (ii„ — i?(w„))^ converges to if g/n ^ and this completes the 
proof of («). 

For proving (ii), recall that 

According to LM2, E{u„) < q'^EjeB^^, (nti(9 " Ij^D) Kj) = 0(^27-^ x 
Therefore, from (3.6) and (3.7), E{\sl\) = E{sl) = E{un) + E{vn) = 
o{'n?^~'^LF'{n)) + 0{q^^~'^L^{q)) and (m) of Lemma 1 follows because q/n van- 
ishes when n goes to infinity. □ 

4. Simulations in dimension d — 2 

The following simulations give an idea of the power of the test under different 
situations of strong dependence. First, one has to approach the asymptotic law 
of the test statistic M„ under the null hypothesis. This is done in subsection 4.1. 
We then focus on the choice of q when n = 128 and n = 256 to guarantee a 
proper size of the test. In subsection 4.3, the power of the test is assessed. 
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Several kinds of long memory random fields are simulated and submitted to the 
testing procedure. As these simulations are highly time consuming, we restrict 
ourselves to random fields of size 128 x 128 and 256 x 256. The simulations results 
reveal a close relation between the power and the kind of strong dependence 
encountered. The power depends on the strength of the long memory but also 
on its anisotropy. 



4.I. Asymptotic law under HO 

The first step to implement the test consists in simulating the asymptotic law 
of M„ under the null hypothesis. According to Proposition 1, this is the law of 

(n'.)^(i)) ^'-[t,.(«<"-(n'-)«<")^' ■ 

(4J) 

where B is the Brownian Sheet on [0, 1]^. After some computations, (4.1) can 
be written 

/ B{ti,t2fdtidt2-2B{1,1) I tit2B{ti,t2)dtidt2 

■/|0,1]2 "'10,1]^ 

-(/ B{tut2)dtidt2] +^^%^/ B{ti,t2)dtidt2 + ^B{1,1)\ 

y\o,ir- J ^ Aoav 

To simulate a sample under this law, each integral above is approximated by a 
Riemann sum, for instance 



ki k2 

n ' n 



[ tit2B{tiM)dtidt2 « ^ V ——B 

where a realization of fi? , . is given by 

sf-,-) =-E E^.-.J., V(fci,fc2)e{l,...,n} 

\ n n I 11 ^ — ^ ^ — ^ 



jl=li2 = l 

with (£j)jgz2 a Gaussian white noise. 

For n = 7000, 10000 reahzations of the law (4.1) have been computed. The 
histogram of the sample is shown in Figure 1. 

From the simulated sample, an estimated mean of 0.0897 is obtained, the 
empirical variance is 0.0018 and the empirical quantilcs of order 90% and 95% 
are respectively 0.1448 and 0.1692. 



4-2. Choice of q 

The sample size n being fixed, one has to choose the value of q involved in 
definition (3.3) of M„. This choice is not easy. We decide to choose it so that 
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0.12 




0.45 



Fig 1. Estimated density function of the limiting law (4.1). 

the size of the test is optimized. This reduces to check that the p-vahies under 
HO are uniformly distributed on [0, 1]. 

SmaU values of q lead to an increase of the probability of rejecting HO. So, 
in order to maximize the power of the test, we have to find the smallest q such 
that the size is correct. 

To achieve this choice, we compute the test for different autoregressive fields 
Xki,k2 defined by 

(1 - aLi)(l - aL2)Xk^^k2 = £kiM^ (4-2) 

where e is a Gaussian white noise, where < a < 1 and Li represents the lag 
operator on the i-th index. The simulation of such an autoregressive field is done 
by filtering a white noise as in [8] . 

Clearly, the size of the test will increase with a. As a consequence, in the 
simulations below, we choose to quote only two situations: a — 0.5 and a = 0.8. 
The last case corresponds to a memory close to strong dependence and a worse 
size in this case might be acceptable. Indeed, it will be impossible to find q such 
that, uniformly on a, the size of the test is strictly the one expected. 

i.2.1. Case n = 128 

Figure 2 represents the empirical distribution of the p- values, computed on 1000 
realizations of A/„, in the case n = 128, when q = 28, q ^ 30 and q = 32. The 
line with crosses stands for a = 0.8 and the line with circles for a = 0.5. The 
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(c) 

Fig 2. Cumulative distribution function of the p-values on [0, 0.1] for model (4-'i) with a = 0.8 
(crosses) and a = 0.5 (circles) when n = 128 and g = 28 (a), g = 30 (b), q = 32 (c). 



diagonal line is added for sake of comparison. The representations are zoomed 
in on [0,0.1]. 

The value (7 = 30 is chosen. For this value, the size associated with a = 0.8 is 
larger than expected (for instance 15% instead of 10%). This error is acceptable 
since the dependence of an autoregressive field (4.2) with a = 0.8 is close to 
long memory. The choice of a larger q reduces this error but on the other hand, 
this may create a bias for smaller values of a (as seen in Figure 2 ior q ~ 32 and 
a = 0.5), and consequently the test becomes less powerful. 

4.2.2. Case n = 256 

Figure 3 represents the empirical distribution of the p-values in the case n = 256 
when q = 35, q = 40 and q = 45. The line with crosses stands for a ~ 0.8 and 
the line with circles for a — 0.5. 

The value g = 40 is chosen for the same reasons as before: The error for 
a = 0.8 seems acceptable in comparison with the loss of power that the choice 
of a larger q would yield. 
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4.3. Power under different alternatives 

We implement the test on different Gaussian long memory random fields. The 
power is assessed according to the type of memory. First, the case when the long 
memory is of tensorial product type is studied: That is when the spectral density 
f{xi,X2) of the field, defined on [— tt, tt]^, is equivalent at to |a;i|"i |2;2|"^, —1 < 
ai < 0, — 1 < a2 < 0. The range of ai and a2 guarantees the integrability of 
/. Notice that this spectral density does not follow the property of Definition 1. 
Then, we report the case when the long memory is isotropic: That is when the 
spectral density is equivalent at to {xf+x'^)"^^ , where, for integrability reasons, 
—2 < a < 0. In the last simulations, we focus on the effect of anisotropy on the 
power. For this purpose, we simulate Gaussian fields whose spectral densities 
are equivalent at to \xi + kx2\°', —1 < a < 0, fc € Z. This is an extreme case 
of anisotropic field since / exhibits only one line of singularity. For all these 
examples, assumption HI holds with 7=1 — (ai + q;2)/2 for the product-type 
case, and 7 = 1 — a/2 for the two other examples (see [14] or [15]). 

All of the above fields are simulated thanks to the spectral method (cf. [13]) 
which consists in the following algorithm: 

1. Generate N independent random variables Zj^'' ),..., (Z^^'', Zj^-*) 

on [— tTjTt]^ according to the spectral measure /i (viewed, up to a normal- 
ization, as a probability distribution); 
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2. Generate N independent random variables Ui, . . . , Un uniformly on [0, 1]; 

3. Compute for all (ij), = ^ J2k=i cos{z[''h + Z^^^ j + 27rL/fc). 

According to the central limit theorem, the resulting field X is Gaussian with 
spectral measure /i when N is large. In practice, the value N — 5000 is fixed. 

4-3.1. Tensorial product type long memory 

The power of the test is assessed for Gaussian fields with a spectral density 
equivalent at to |a;i|"i|x2|"^, — 1 < ai < 0, — 1 < a2 < 0. Figure 4 shows a 
simulation of such fields on a 256 x 256 grid using the spectral method, when, 
from left to right, ai = a2 = —0.25, ai = 02 = —0.5 and ai = a2 = —0.75. 
These cases correspond respectively in (2.5) to 7 = 1.25, 7 = 1.5 and 7 = 1.75. 

The empirical c.d.f. of the p-values of the test is represented on Figure 5. 
Each curve is computed on 500 simulated fields. On the left, the simulations 




Fig 4. Gaussian fields with product-type long memory where 7 = 1.25 (left), 7 = 1.5 (middle) 
and 7 = 1.75 (right). 




Fig 5. C.d.f. of the p-values when 7 = 1.25 (crosses), 7 = 1.5 (circles) and 7 = 1.75 (trian- 
gles), where n = 128 (left) and n = 256 (right). 
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correspond to fields of size 128 x 128, wlrile on tlie riglit n — 256. The parameter 
q has been chosen according to subsection 4.2, that is g = 30 when n = 128 and 
q = AO when n = 256. In each case, three curves are represented, corresponding 
to ai — a2 = —0.25 (crosses), ai = a2 = —0.5 (circles) and ai = a2 ~ —0.75 
(triangles). We observe logically that the power of the test increases with the 
strength of the memory (quantified by 7) and with the size of the sample. 

4-3.2. Isotropic long memory 

Now, the power of the test is assessed on isotropic long memory fields. Their 
spectral density is equivalent at to {x^ + y^)"^^ , —2 < a < 0. A simulation 
of such fields is presented on Figure 6 when a = —0.5, a = —1 and a = —1.5, 
which correspond to 7 = 1.25, 7 = 1.5 and 7 = 1.75 in (2.5). 

Figure 7 represents the same c.d.f as in Figure 5 but for the isotropic case 
when a = —0.5 (crosses), a = — 1 (circles) and a = —1.5 (triangles). These 




Fig 6. Gaussian fields with isotropic long memory when 7 = 1.25 (left), 7 = 1.5 (middle) 
and 7 = 1.75 (right). 




Fig 7. C.d.f. of the p-values when 7 = 1.25 (crosses), 7 = 1.5 (circles) and 7 = 1.75 (trian- 
gles), where n = 128 (right) and n = 256 (left). 
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choices correspond to the same strengths of long memory (i.e. the same 7) as in 
Figure 5. The power follows a similar behavior: It increases with the strength 
of the memory and with the sample size. However, the power for the isotropic 
case is smaller than in the product-type setting studied before. This is due to 
the form of the statistic Af„ (see Section 3.1): The empirical variance of 5 * ^ 
is computed using quadrants Aj^s. This suits better product-type fields than 
isotropic ones for detecting long memory. Indeed, in the product-type setting, 
this empirical variance will tend to return higher values. This sensitivity to 
anisotropy is studied further in the last subsection. 

4-3.3. The effect of anisotropy on the power 

As seen before, the power for the isotropic long memory case is smaller than 
for the product-type setting. To assess properly the sensitivity to anisotropy, we 
focus on one-direction long memory fields in the sense that their spectral density 
behaves at zero as \xi +kx2\°', k £ Z. Note that this form of the spectral density 
demands —1 < a < to be integrable. This yields the restriction 7 < 1.5 in 
(2.5). Figure 8 shows a simulation on a 256 x 256 grid when fc = — 1 and fc = 
and when a = —0.5, corresponding to 7 = 1.25 in (2.5). As we can see on 
Figure 8, the worst case, in terms of computation using quadrants, should be 
fc = — 1 while the most suitable situation should be fc = 0. 

Figure 9 shows the power of the test when k = ~1 (solid line) and fc = 
(dotted line) for n = 128 (left) and n — 256 (right). This representation confirms 
the effect of anisotropy. Therefore, from a practical point of view, it seems better 
to study first the isotropy of the dependence before testing the presence of long 
memory, in order to rotate the image sample if necessary. 




Fig 8. Gaussian fields with one-direction long memory (7 = 1.25) when k = —1 (left) and 
k = Q (right). 
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Fig 9. C.d.f. of the p-values when k = —1 (solid line) and k = (dotted line) for 7 = 1.25 
and where n = 128 (left) and n = 256 (right). 
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